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It has recently been emphasized that all known exact evaluations of gap probabilities for 
classical unitary matrix ensembles are in fact r-functions for certain Painleve systems. We 
show that all exact evaluations of gap probabilities for classical orthogonal matrix ensembles, 
either known or derivable from the existing literature, are likewise r-functions for certain 
Painleve systems. In the case of symplectic matrix ensembles all exact evaluations, either 
known or derivable from the existing literature, are identified as the mean of two r-functions, 
both of which correspond to Hamiltonians satisfying the same differential equation, differing 
only in the boundary condition. Furthermore the product of these two r-functions gives the 
gap probability in the corresponding unitary symmetry case, while one of those r-functions 
is the gap probability in the corresponding orthogonal symmetry case. 



1 Introduction 

An ensemble of N x N random matrices X with joint probability density of the matrix elements propor- 
tional to 

oo N 

cxp(5> 3 Tr(X->)) -Hgixj), (1.1) 

Xj denoting the eigenvalues, is invariant under similarity transforms X i— > A~ 1 XA. In particular, if X 
is an Hermitian matrix with real, complex and quaternion real elements, labelled by the parameter /3 
taking the values p = 1, 2 and 4 respectively, then the subgroups of unitary matrices which conserve this 
feature of X under similarity transformations are the orthogonal (/3 = 1), unitary (/3 = 2) and unitary 
symplectic matrices (f3 = 4). For this reason the ensemble is said to have an orthogonal (j3 = 1), unitary 
(j3 = 2) or symplectic symmetry ((3 — 4). The eigenvalue probability density function (PDF) for these 
ensembles has the explicit form 

^n^) n i^-^r, (i-2) 

j=l l<j<k<N 



The function g(x) in (LI) and (L2) is referred to as a weight function. In the cases j3 = 1 and [3 — 2 
the weight functions 

e -x 2 /2^ Gaussian 
, _ . a>-i)/2 e -*/2 ( x > o), Laguerre 
gW- gi (x)-l (1 _ x)(o _ 1)/2(1 + x) {b-i)/2 (_ 1<a . < i) j j acobi V- 6 > 



(1 + .t 2 )-( q+1 )/ 2 , Cauchy 



and 



e x , Gaussian 
, J xa e~ x {x>0), Laguerre 

9^=^ = < (1 -*).(! + (-!<*<!), Jacob! ,L4) 
[l + x 2 )~ a , Cauchy 



are said to define classical matrix ensembles with an orthogonal and unitary symmetry respectively, or 
simply classical orthogonal and unitary ensembles (a similar definition applies in the symplectic case — 
see e.g. 0). We recall (see e.g. the introduction of jjJJ) that the Cauchy ensemble includes as a special 
case the PDF 

1 -Q | e *WY ; (L5) 

l<j<k<N 

which specifies the Dyson circular ensembles. The two ensembles are related by the stereographic pro- 
jection 

= jj— (1.6) 



In particular, changing variables in (1.5) according to (1.6) gives a PDF of the form (1.2) with g(x) a 
Cauchy weight function, which in the cases (3 = 1 and 2 is specified by ([ll]) and (|0|) with a = N. 

Our interest is in a special property of the probability Ef}(0; I; g(x); N) of having no eigenvalues in the 
interval / when the eigenvalue PDF is specified by fll-2| ) in the case that g (x) is classical. The probability 
is specified as a multiple integral by 

1 N f 

E (O;I;g(x);N) = -l[ dx 3 g( Xj ) ]J \x k - Xj f , (1.7) 

J = 1 JlaV l<j<k<N 

where Iq is the interval of support of g(x). The special property is that for g{x) classical Ep admits 
Painleve transcendent evaluations for certain / (the evaluations are in some cases restricted also to certain 
scaled limits). We will focus on a structural aspect of these formulas, by showing that in the orthogo- 
nal case all known Painleve transcendent evaluations can be identified as r-functions for Hamiltonians 
associated with the Painleve functions, and in the symplectic case as the mean of two r-functions. 

Our work builds on the recently emphasized |l3|, Q fact that all gap probabilities for classical unitary 
ensembles that have been characterized as the solution of a single differential equation, are in fact r- 
functions for certain Painleve systems. Such characterizations of the gap probability for classical unitary 
matrix ensembles are known when the gap consists of a single interval including an end-point of the 
support |27| |2[ [l7| |l[ (7) , or a double interval symmetrically placed about the origin again including the 
end-points of the support or the origin (applicable to even weight functions only) |27[ |3^, pit - In the 
special case of the gap probability for scaled, infinite GUE matrices in the bulk, the identification as a 
r-function for a Painleve system was made by Okamoto and quoted in the original paper of Jimbo et 
al. |2^, pg. 152] deriving the Painleve evaluation. For the more general problem of characterizing the 
gap probabilities in the case of multiple excluded intervals, the fact that the probability is the T-function 
for certain integrable systems associated with monodromy preserving deformations of linear differential 
equations with rational coefficients was a main theme of p(j) , and then generalized to a more general 
setting (but not the most general case of interest in random matrix theory) by Palmer p5| . Harnad 
and Its |l8[| have recently discussed the work of Palmer from a Riemann-Hilbert problem perspective. 
Identifications of the gap probabilities in the case of multiple excluded intervals as r-functions in the Sato 
theory is a theme of the work of Adler, van Moerbeke and collaborators (see e.g. ||). 

The situation with the exact evaluation of gap probabilities for matrix ensembles with an orthogonal 
symmetry is immediately different due to the restricted number of evaluations in terms of Painleve 



transcendents presently known [gCj [10 , 1 1 1 . In the orthogonal case, the exact evaluations can be catalogued 
into two distinct mathematical structures — the finite N ensembles and their scaling limit for which the 
r-function identification is immediate, and the infinite Gaussian and Laguerre ensembles scaled at the soft 
and hard edges respectively in which the known Painleve transcendent evaluations reduce to a r-function 
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92 (x) 


x- 1 ' 2 g 2 (x 1 l 2 ) 


classical ensemble 


-x 2 

e x 
(l-x 2 ) a 
(l + x 2 )- a 


x ~l/2 e -x 

x- l / 2 {l - x) a 


Laguerre: a = — 1/2 
Jacobi: x \— > 2x — 1, b = — ^ 
Jacobi: x i-> b = -1/2, a = a - N + \ 



Table 1: Even classical weights, their transformed form, and the corresponding classical unitary 
ensemble with N/2 eigenvalues. 

after some calculation. In the symplectic case all known exact evaluations result from a formula relating 
the gap probability in the symplectic case to that in the orthogonal and unitary cases. Further special 
features of the exact evaluations in the orthogonal and unitary cases then allows the exact evaluations 
in the symplectic case to be identified as the mean of two r-functions, both of which correspond to 
Hamiltonians satisfying the same differential equation, differing only in the boundary condition. 



2 Orthogonal matrix ensembles 
2.1 Finite N ensembles 

It has been shown in |1(J that for the classical weights (1.3), having the additional property of being even 
(which is the case for the Gaussian, symmetric Jacobi (a = b) and Cauchy weights), 



(2.1) 



i/ 2)i 



E 1 (0; (s, s); gi (x);N) = E 2 (0; (0, s 2 ); x' 1 ' 2 g 2 {x 1 / 2 )- N/2) 

where on the RHS x > 0, and it is assumed N is even. Now a unitary ensemble with weig htx- 1 / 2 g 2 [x 
in which g 2 (x) is an even classical weight, is equal to another unitary ensemble with a classical weight, 
after a suitable change of variables as detailed in Table [jj. Hence it follows that 

E 2 (0;(Q lS 2 );x^ 2 g 2 (x^ 2 );N/2) 

E 2 (0;{0,s 2 );x- 1/2 e- x ;N/2), Gaussian 
# 2 (0;(-l,2s 2 - l);(l + a;)- 1 ' /2 (l-a;) a ;iV/2), symmetric Jacobi (2.2) 

^ E 2 (0- (-1, (s 2 - l)/(s 2 + 1)); (1 + x)- 1 ' 2 ^ - x) a - N+1 / 2 ;N/2), Cauchy 

This substituted in ( |2.l|) gives -Ei(0; (— s, s); gi(x); N) for the even classical orthogonal ensembles in terms 
of E 2 for certain classical unitary ensembles. The latter furthermore have the gap free interval including 
an end-point of the support of the weight function. In such a case, we can deduce from the existing 
literature that E 2 , and consequently E\, is a r-function for an appropriate Painleve system. 

Consider first ^(O; (0, s); x a e~ x \ N), specifying the probability that there are no eigenvalues in the 
interval (0, s) of the Laguerre unitary ensemble. Following [^3] and 1 19 introduce the Hamiltonian Hy 
associated with the Painleve V equation by 



tHy = q{q - l) 2 p 2 - {(«a - Vi)(q ~ l) 2 - 2(«i + v 2 )q(q - 
where the parameters v\ , . . . , are constrained by 

V\ + v 2 + v 3 + v 4 



1) + tq}p + (w a - vx)(v 4 - vi)(q - 1) (2.3) 



0. 

The relationship of ( |2.3D to Py can be seen by eliminating p in the Hamilton equations 

, _ aff , dH 

9 dp ' 



P 



(2.4) 



(2.5) 
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One finds that q satisfies the equation 



y = [kt. 



1 



2y y-i 



(//' 



(y-i) 



a V H — ) + 7t + ^ 



with 



1 



i 



"=0(^3-^4), /3 = -r(«2-«i) , 7 = 2vi+2v 2 -l, 5 = — — ■ 



(2.6) 



(2.7) 



This is the general Py equation with 6 = — i (recall that the general Py equation with 5^0 can be 
reduced to the case with S = — \ by the mapping t > \J — 2St). Now introduce the auxiliary Hamiltonian 



(T V = tHy + (l>3 - Vt)(V4 - Vi). 



(2.8) 



Of course with tHy replaced by cry in (2J)), the Hamilton equations remain unchanged so ay /t is also 
a Hamiltonian for the same Py system. The quantity cry satisfies the second order, second degree 
differential equation 

(ta"f -[a- to' + 2(a') 2 + K + v\ + v 2 + ^K] 2 + 4(^ + <r')Oi + a')(v 2 + a')(v 3 + a') = 0, (2.9) 
with 

^o = 0, vi=v 2 -vx, v 2 = v 3 -v 1 , v s =v 4 -Vi (2.10) 

(because (|2.9|) is symmetrical in {vk} any permutation of these values is also valid). Conversely, each 
solution (with a" 7^ 0) of (2.9) leads to a solution of the system (2.5) |23fl . 
The r-function associated with the Hamiltonian ay /t is specified by 



ay = t—\ogT av (t). 
But from the work of Tracy and Widom [E7| we know that 



t-logE 2 (0;(0,t);x a e-*;N) 



(2.11) 



(2.12) 



satisfies (2.9) with 

vq = 0, ^i = 0, v 2 — N + a, 1/3 = N, 
subject to the boundary condition 

a{t) 



T(N + a+l) ta+ ^ 



t^o+ r(7V)r(a + l)r(a + 2) 



Consequently, after equating ( 2.11 ) and ( 2.12 ), and normalizing r CTv so that t<7 V (0) = 1, we have 



E 2 (Q-{Q,t)-x a e~ x -N) =r <Ty (t) 



"0=0. "1=° 
u 2 = N + a, u 3 =N 



(2.13) 



With a = —1/2 we see that this corresponds to the Gaussian case of (2.2). Recalling (2.1) then gives the 
sought r-function formula for the Gaussian orthogonal ensemble, 



E 1 (0;(-s,s);e- x2 / 2 ;N)=r av (s 2 ) 



"0 = 0, "1=0 
" 2 =JV-l/2, " 3 = ]V 



(2.14) 



Consider next E 2 (0: (— 1, s); (1— x) a (l+x) b ; N), specifying the probability that there are no eigenvalues 
in the interval (— l,s) of the Jacobi unitary ensemble. According to (2.2) this is relevant to both the 
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symmetric Jacobi and Cauchy cases. Introduce the Hamiltonian ifyi associated with the Painleve VI 
equation by p3| 

t(t - l)H vl = q(q -l)(q- t)p 2 - { X o(q - l)(q - t) + xMQ ~ t) + {6 - l)q(q - l)}p + x(Q - t) (2.15) 
where 

x = \{xo + xi + o-i? - \xlo- 



Eliminating p from the corresponding Hamilton equations (2.5) shows that q satisfies the Pyi equation 

„ i n 



y =o - + — 

2 Vy y 



L_ + JLW)* _ (i + JL + JL_V 

-1 y-t) yy> \t t-1 y-tJ y 



y(y- i){y-t) 
t 2 (t-i) 2 



t t-i J(t-i) 



with 



Furthermore, the auxiliary Hamiltonian 

hvi = t(t - 1)H VI + (6163 + hh + b 5 bi)t - - b i b ^ 

where 



(2.16) 



h = ^{xo+Xi), b 2 = -{Xo-Xi), & 3 = ^( 



l<j<fc<4 
1 



1 + Xoo), h = -(9 - 1 - Xoo), 



satisfies the differential equation 

4 

h! vl (t(l - i)/i VI ) 2 + {h' Yl [2h m - (2t - l)/i VI ] + hb 2 b 3 b 4 } 2 = JJ(ftVi + b 2 k ) 



(2.17) 



k=l 



and conversely, each solution of (2.17) such that /i VI 7^ leads to a solution of the corresponding Hamilton 
equations. Now, we know from the work of Haine and Semengue ]l7| ], and Borodin and Deift JtJ, that 



a(t) := t(t - 1)1 log£ 2 (0; (-1, -1 + 2t); (1 - x) a {\ + x) b ;N) - b x b 2 t + ^(hb 2 + 6364) 



with 



61 = 6 2 = N 



a + b 



b* = 



a + b 



b A = 



- b 



satisfies ( 2.17 ). Comparing with ( 2.16 ) we see that with this choice of parameters 

t{t - l)llog£ 2 (0; (-1, -1 + 2t); (1 - x) a (l + x) b ;N) = h vl + b x b 2 t - ^{b x b 2 + b 3 b±). (2.18) 



Thus, denoting the RHS of ( 2.18 ) by hvi, we see from ( 2.16 ) that hyi/t(t — 1) is a Hamiltonian for the 
Pyi system, and defining the corresponding r-function by 



we have that 



hvi=t(t-l)-\ogT~ hvi (t) 



E 2 (0; (-1, -1 + 2t); (1 - x) a (l + x) b - N) = rtti) 



b 1 =b 2 =N + (a + b)/2 
b 3 = (a + b)/2, b 4 = ta-b)/2 



(2.19) 



Recalling ( pT2| ) and (2.1) then gives the sought r-function formulas for the gap probabilities in the Jacobi 
orthogonal and Cauchy orthogonal ensembles, 

E x (0;(-s,s);(l- X 2 )^ 2 ;N) =r~ h (s 2 ) 



E x (0; (s, s); (1 + z 2 )-^ 1 )/ 2 ; N) =r^ (-^) 



i. 1 = i. 2 =-N/2+(<l-l/2)/2 
b 3 =(a-l/2)/2, 6 4 = (a+l/2)/2 

2 



b 1 =b 2 = a/2 
6 3 = (o-JV)/2, 6 4 = (a-JV + l)/2 



(2.20) 
(2.21) 
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2.2 Bulk scaling limit 

Let us consider now the N — ► oo bulk scaling limit of an orthogonal ensemble, and the quantity 
£ , J )ulk (0; 2s) specifying the probability that there are no eigenvalues in an interval of length 2s with 
the mean spacing between eigenvalues equal to unity. By an appropriate scaling, each of the probabilities 
in ( 2.14 ), ( 2.20 ) and ( 2.21 ) tends to £^ ulk (0; 2s). For example, in the Gaussian case the required scaling 
is s i — ► 7ts/V2-/V and so 



^ bulk (0; 2s) = lim EAQ-A 



);e' x ' 2 ;N). 



This scaling applied to (2.14) is known to lead to the result [jlO| 

^ ulk (0;2s) = exp(-^ 7r " a B (t) 



dt 

a=-l/2 t 



where <7E?(t) satisfies the equation 



(ta' B ) 2 + a' B (a B - ta' B )(4a' B - 1) - a 2 {a' B f = 



(2.22) 



(2.23) 



subject to the boundary condition 

1 



J 2 (Vt)~ J a +l{Vi)Ja-l(y/i) 



t^o+ 2 2 + 2a r(l + a)r(2 + a)' 



(2.24) 



In fact the expression ( 2.22 ) is precisely the r-function for a particular Pni system. To see this, 
following Okamoto introduce the Hamiltonian 

tH = q 2 p 2 - (q 2 + viq - t)p + ^-(vi + v 2 )q. 



Substituting this form of H in the Hamilton equations ( |2.5| ) and eliminating p shows that y(s) = q(t)/s, 
t = s 2 , satisfies the general Painleve III equation (Painlcve III' in the notation of p4j) 



d 2 y = }_/^y\ 2 _ l^y + lr ay z + p\ 

ds 2 v\ds J s ds s 



'72/ + - 
V 



with 



a = — 4u2, /3 = A(vi + 1), 7 = 4, S = -4. 
It is shown in |24j| that the auxiliary quantity 



satisfies the equation 



1 , 1 

h = tH + -v? * 

4 1 2 



(th") 2 + Vl v 2 ti - (4(ti) 2 - l)(h - th') - + v 2 ) = 0, 



(2.25) 



(2.26) 



and conversely all solutions of this equation (assuming h" ^ 0) lead to the Pnr system. It is a simple 
exercise to verify from the fact that h satisfies ( 2.26| ), the result that 



<rm>(fi) := -(tH) 



- -r( v i ~ v 2) + -; 
4 4 



(2.27) 



satisfies the equation 



Kir) 2 - «iwa(om0 2 + °m>(4*>m> - l)(^ni' - ta' liv ) - ^ («i - v 2 f = 



(2.28) 
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We note from (2.27) that —atu'{t)/t is a Hamiltonian for the Pur system, so we can introduce the 
corresponding r-function by 



aw (i) 



-t^logmHt). 



(2.29) 



Now the equation (2.28) with v\ = v-i = a is identical to ( [2.23] ), so comparison of ( |2~29| ) and $L2% ) gives 
the r-function evaluation 



i? 1 bulk (0;2 S ) = r II p(7r 2 ^. 

v l=v?=— 1/2 

The boundary condition satisfied by cthf(£) is the a = — 1/2 case of ( |j4[ ), 



CTiii'(i) ~ TT - 

t^0+ 27T 



3in2Vt' 



2y/i 



t— >0+ 7T 



(2.30) 



(2.31) 



2.3 Cumulative distribution of the largest eigenvalue in the scaled infinite 
GOE 

The GOE has the property that to leading order the support of the spectrum is confined to the interval 
[—y/2N, \/2N]. It was shown in || that by scaling the eigenvalues 



Ai 



'2N - 



X 



s/2NV 6 ' 



(2.32) 



so that the origin is at the right hand edge of the leading support and the eigenvalue positions then 
measured in units of 1/ V2-/V 1 / 6 , the distribution functions describing the eigenvalues in the neighbourhood 
of this edge (referred to as a soft edge since the density on both sides is non-zero) are well defined. 
It was shown by Tracy and Widom [Sol (see [O for a simplified derivation) that 



i?r ft (0; (a, oo)) := F^s) := lim eJo; 



e 2 



'2N ■ 



e 2 



V2NV e 



,oo);N 



where q(t) is the solution of the non-linear equation 



tq + 2q i , 



subject to the boundary condition 



q(t) ~ -Ai(t) 



as 



t — > oo, 



(2.33) 



(2.34) 



(2.35) 



where Ai(t) denotes the Airy function. (Here we have replaced q by —q relative to its use in the original 
work; this is valid because ( 2.34 ) is unchanged by this mapping.) Since the general Painleve II equation 
reads 



tq + 2q A 



a. 



(2.36) 



( p.34[ ) is the special case a = of Pn • Thus ( 2.33 ) represents an explicit evaluation of the gap probability 
in terms of a Painleve transcendent. It is the objective of this subsection to show that in fact ( p. 35 ) can 
be identified as a r-function corresponding to the Painleve II system with a = 0. Consequently its 
logarithmic derivative satisfies a single nonlinear differential equation. 

Now, in the case of the probability analogous to F\ (s) in the infinite, scaled Gaussian unitary ensemble 
(GUE), the known exact evaluation pq| allows one to immediately make an identification with a r-function 
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|T3[ . It is relevant for the purpose of identifying (2.33) to revise the theory underlying this result. Tracy 
and Widom |2£| have derived the result 

£| oft (0; (s,oo)) := F 2 (s) := lim E 2 (0; {V2N + s/V2N 1/6 , oo); AT) = exp f - / R{t)dt), (2.37) 

where satisfies the second order second degree differential equation 

(R") 2 + AR'((R') 2 - tR' + R) = 0, (2.38) 

and have furthermore derived the alternative formula 

F 3 {a) = e-^*-'^ dt ) (2.39) 

where q(t) is the same Painleve II transcendent as in ( [2.33 ). 

To see how the evaluation ( 2.37 ) relates to a r-function for the Painleve II system, we recall that in 
the Hamiltonian formalism of the Pn equation [Ell , one defines a Hamiltonian Hjj by 



H u = --(2q 2 - P + t)p-( a + -)q. (2.40) 

The canonical coordinate q and momenta p must satisfy the Hamilton equations (|2.5| ) . Elimination of the 
variable p between these equations shows that q satisfies the Painleve II equation (2.36). Furthermore 
the Hamiltonian ( 2.40 ), regarded as a function of t, satisfies the second order second degree differential 
equation 



(i^) 2 + 4(ff XI ) 3 + 2H' n {tH' n - Hu) - \{a + \f = 0, 



(2.41) 



referred to as the Jimbo-Miwa-Okamoto a form for Pn . It is also straightforward to show that H can 
be expressed in terms of the Painleve II transcendent q according to 



Hu = \{q'?-\^ + \t) 2 -{a + \)q. 
Finally, we recall that the r-function associated with the Painleve II Hamiltonian is defined by 



(2.42) 



-fr" = 37 log ''"li- 
ra 



Setting 



u(t;a + -) 



-2 1/3 H u {-2 1/3 t) 



we see from ( 2.41 ) that u satisfies the equation 



(2.43) 



(2.44) 



(u") 2 + 4m' [(u') 2 - tu' + u] - (a + -) 2 = 0. 
Comparison of ( 2.45| ) with (2.38) shows 

R(t) =u(t;0) = -2 1/3 ff n (-2 1/3 i)| _ i a< 



In light of this identification, comparison of (2.37) and (2.43) then shows, 

F 2 ( S )=r„(-2 1 / 3 S ) 



= -1/2 



(2.45) 



(2.46) 



(2.47) 



The appropriate boundary condition for this r-function is most simply expressed in terms of R(t), 

R(t) ~ [Ai'(t)] 2 -t[Ai(t)] 2 . (2.48) 

t — >OQ 
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Thus, up to a scale factor, -F^s) is precisely the r-function associated with the Hamiltonian (2.40) for 
the Painleve II system with a — —1/2. A curious feature of (2.46), which follows from (2.42), is that 
R(t) is naturally expressed in terms of the Painleve II transcendent q = q(t; —1/2), whereas the result 
( [2.39D involves the Painleve II transcendent with a = 0. In particular, (2.37) and ( p. 39 ) give 



R'(t) = -q 2 (t;()) 

while (2.46), ( 2.40| ) and the first of the Hamilton equations (2.5) give 

R 'V = "2^3 -V2) + Q 2 (t, -1/2) + t 

In fact, as noted in [|13|] , for e = ±1, it is true that jl6) 

_ e2 V3 g2( _ 2 -i/a t) o) = | ?(t) l e) _ eq 2 {ti l e) _ 1 et) 



tl ^_2 1 /3 t 



(2.49) 



(2.50) 



(2.51) 



which reconciles (2.50) with (2.4E). 

We are now in a position to identify (2.33) with a r-function. The formula Q2.50 ) is just the special 
case a — of the identity 



dt 



H u (t) 



ct=a-l/2 



dt v ' ; 



q'(t,a- 1/2) + q 2 (t,a- 1/2) + 



(2.52) 



which is derived from (2.44), ( [2.4C ) and the first of the Hamilton equations (2.5). To make use of this 

(2.53) 



result we first note that the equation (2.33) can be written 

Fi(s) = e - 5L 00 ( t -*X? 2 (*)+s'(t))<*t 



The identity (2.52) with a — 1/2 allows this in turn to be rewritten as 



Fi (s) = exp 



> (t- S )|[2- 1 /3 u( _2-V3 i;1/2) _^ 



exp 
exp 



2- 1 / 3 U (-2- 1/3 t;l/2)-^ 
8 



dt 



H n (t) 



a=0 8 



dt 



(2.54) 



where the final equality follows from ( 2.44 ). 

We now associate with H the auxiliary Hamiltonian 

t 2 

h u = H n + -. 



(2.55) 



Of course, the Hamilton equations (2.5) remain valid for H replaced by h, so h is also a Hamiltonian for 
the same Painleve II system. Introducing the corresponding r-function by 



hi = ^ l °S T hn, 



we see from fl2.56|) that 



Fi{s) = T hn {s) 



a=0 



(2.56) 



(2.57) 



which is our sought result. Note that h\\ satisfies ( [2.4l| ) with the substitution Hu = h\\ — i 2 /8. It follows 
from ( 2.42| ), ( |2.55 ) and (2.35) that we seek the solution of this equation with a — and such that 



hu(t) ^ ^Ai(i) + \ {[M'(t)] 2 - t[M(t)] 2 } 



(2.58) 
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Unlike the situation with Sf ulk (0;2s), in which the corresponding finite system gap probability 
£1(0; (— s, s); e~ x ^ 2 ;N) is itself a r-function, there is no known Painleve transcendent evaluation of 
the finite N quantity in the definition ( 2.33| ) of Fi(s). Nonetheless, ( 2.57 ) can be obtained as a limiting 



sequence of finite N Painleve transcendent evaluations, which in fact is how we were led to (2.57) in the 
first place H]. The finite N results are not for gap probabilities though^. Rather they relate to the 
quantity fm specifying the number of fixed point free involutions of {1, 2, ... , 2./V} constrained so that 
the length of the maximum decreasing subsequence is less than or equal to 21. This is specified by the 
generating function 

which from the work of Rains |2(| (see also ](|) has the integral representation 
-i l /2 , i , l fit fir t 1 

j=i i<j<k<l 

where Zj :— e l9j . Although not at all obvious from the definition, it has been proved in || that 

limP^O- ^/2) 1/3 )) =Fi(s). (2.61) 

I— >oo Z 

The significance of this result from the present perspective is that we have recently shown Jl4|] Pi (t) to be 



p ^=^r-(-k) 1 f r de ie 2t ^^f[\i- Zj \ 2 n \i- Zl z k \ 2 \z 3 ~ Zk \ 2 : (2.60) 



equal to the r-function for a certain Painleve V system which scales to the result (2.57) (the evaluation 



of Pi{t) in terms of a transcendent related to Painleve V was first given by Adler and van Moerbeke Q). 

2.4 Cumulative distribution of the smallest eigenvalue in the scaled infinite 
LOE 

In the LOE, as N — > oo the spacing between the eigenvalues in the neighbourhood of the origin (referred 
to as the hard edge because the eigenvalue density is strictly zero for x < 0) is of order 1/N. With the 
scaling 

the distribution functions describing the eigenvalues near the hard edge have well defined limits ||. Our 
interest is in 

< rd (0;(0, S );(a-l)/2):= lim ^(0; (0, a^" 1 )/^-*/*; N ), 

which is equal to the probability of no eigenvalues in the interval (0, s) of the scaled, infinite LOE, or 
equivalently to the cumulative distribution of the smallest eigenvalue in the ensemble. It has been shown 
to have the Painleve transcendent evaluation []ll| 



£f rd (0; (0, s); (a - l)/2) = exp ( - ± J' ( log |)« 2 (t) dt) exp ( - ± f ^ dt) , 
where q(t) satisfies the nonlinear equation 



(2.62) 



t{q 2 - l)(tqj = q{tq'f + \{t - a 2 )q + ^tq 3 (q 2 - 2). (2.63) 



1 Since completing this work the gap probability i?i(0; (s, oo); e x ^ 2 \ N) has been evaluated as a Pv r-function 
|15|. and it scales to Fi(s). 
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This equation, which is to be solved subject to the boundary condition 

1 



a/2 



is transformed |29f| via the substitutions 



t = x 2 , q(t) = 
to the Py equation ( |2.6| ) for y(x) with parameters 



1 + y{x) 
1 - y{x) 



-2. 



(2.64) 



(2.65) 



(2.66) 



In this subsection we will show that (2.62) can be identified with a r-function corresponding to the 
Py Hamiltonian (2.3). 

To begin we observe that 



in which use is made of (2.64) for its derivation. Hence we can write 

1 r 

8io 



(logs - logt) 



^ ard (0;(0, S );(a-I)/2) = exp( 

/ 1 f vs 
= 6XP V~ 4 J ( l °Ss-^ogt) 



t- 1/2 q + 2t 1/ ~q' 



2 ' 



dt 



dq , 
x— + q + xq 
dx 



dx 



(2.67) 



But it follows from (2.65) that 



dq 2 dy 

dx (1 — y) 2 dx ' 



and thus 



x— + + x 2 = 1 
dx (1 — y) 2 dx 



(2a: ^ - y 2 + Axy + 1) + x. 



Consider now the Hamiltonian (2.3). With the replacements 



1 ^ 



q h- > y, pHZ, t 1 ► TjX, Hy 1 ► — i?y 



(2.68) 



(2.69) 



it reads 



xH v (y,z) = y{y - l) 2 z 2 - |(t) 2 - wi)(y - l) 2 - 2(«i + ^2)2/(2/ - 1) + nxy\z 

+ (v 3 -v 1 )(vA-vi)(y-l) (2.70) 

According to (2/7), the remark below (2.7) in parenthesis and (2^4), the parameter values (|2.66| ) correspond 
to the Hamiltonian (2.7C) with 



?7 = 2, vi = -v 3 = -7(0- l)j «2 = -U4 = \{a + 1) 



(2.71) 



Furthermore we need to add a term — l(a 2 — 1) to the Hamiltonian in order that a well-defined limit for 
the auxiliary Hamiltonian specified below exists as x — > + . Making use of the Hamilton equations it 
follows that with these parameter values 

d 



dx 



(xHy) — —2yz 

x ix~ = 2y(3) ~ 1)2 z ~ i (y ~ 1)2 + v(y ~ x) ~ 2xy - 



(2.72) 
(2.73) 
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Substituting for yz in (2.73) using (2.72) we see that 

1 L dy 9 

2x-f- - y + Axy + 1 

dx 



(i - y) 2 



= -2 



dx 



(xHv) 



a - 1 



Substituting this in (2.68) then gives 

dq 

x—+q^ 

dx 



xq 2 = -2 



a - 1 1 



dx 



Finally, substituting ( 2.74 ) in ( 2.67 ) and integrating by parts we arrive at the result 



E 



hard 



(0;(0, S );(a-l)/2)) = exp / ( 



a- 1 1 



- x + Hy ) dx 



2 4 

Thus if we define the auxiliary Hamiltonian and corresponding r-function for the Py system by 



r t-t 1 o-l 
hy = Hy - -x H — . 



kv= dx- l0ST ^ 



we obtain the sought r-function evaluation 

^ ard (0;(0, S );(a-l)/2)=r, (y^) 



V =2 

i>i=-D 3 =— -j(a-l) 
■U2=-«4 = j(a+l) 



Note that with the parameters (2.71) it follows from (2.8) and (2.69) that 

xhy 



1 2 (a - 1) a(a- 1) 
M 1 )-^ + — 2 — 4 : 



where crv(x) satisfies (2.9) with 1 1— > 2x. The boundary condition for this Hamiltonian is 



x -o^W ~ 7 x2 [ J a( x ) " J +i(a;)Ja-i(a;)] ■ 
x^o+ 2 4 



(2.74) 



(2.75) 



(2.76) 



(2.77) 



(2.78) 



(2.79) 



The r-function evaluation of £^ ard differs from those of E^ and Ef in that no finite N quantity 
is known which is itself a r-function, has an interpretation as a probability, and which scales to (2.77). 



3 Symplectic matrix ensembles 
3.1 Finite TV ensembles 

With TV finite there is in fact only one symplectic matrix ensemble — the circular symplectic ensemble — 
for which the gap probability can be written in terms of Painleve transcendents using results known in the 
literature^]. With EpiO; (—(f>, 4>); N) denoting the probability that there are no eigenvalues in an interval 



(— <ft, (f) of the circular ensemble specified by the PDF (1.5), this is possible due to the inter-relationships 
between gap probabilities due to Dyson and Mehta |8| 

implying the evaluation of E4 from knowledge of the evaluation of E\ and E 2 



Regarding the latter, let 4> be related to s via the stereographic projection formula (1.6) with 9 
x s. Then from the relationship between the circular ensemble and Cauchy ensemble we have 

E 1 {Q; (-0, 0); 27V) = ^(0; (-*, s); (1 + x 2 )^ N+1 ^ ;2N) 
E 2 (0; (-</), <j>); 27V) = E 2 (0; (s, s); (1 + x 2 y 2N ; 27V). 



2 Since completing this work the gap probability E4(0; (s,oo);e s ; N) has been evaluated as the arithmetic 
mean of two Pv r-mnctions, the Hamiltonians of which satisfy the same differential equation. 
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Now we know from (2.1) that 

^(0; (s, s);(l + x 2 )- (2N+1)/2 ; W) = E 2 (0; (0, S 2 ); x^l + x)- 2N ;N) 
while an identity in |]lof gives 

E 2 (0;(-s,s);(1 + x 2 )- 2N ;2N) 

= E 2 (0; (0, s 2 ); x~ l l 2 (l + x)~ 2N ■ N)E 2 (0; (0, s 2 ); x x l 2 (\ + x)~ 2N ;N). (3.2) 

Thus we have 

E 4 (0; (-0, <f>);N) = \{e 2 (0: (0, s 2 ); ^(l + x)~ 2N ;N) + E 2 (0; (0, s 2 ); x^ 2 (l + x)~ 2N ; N)) 
= \{e 2 (Q: (-1, (s 2 l)/(s 2 + 1)); (1 + x)^' 2 (l xf/ 2 ;N) 
+ E 2 (Q: (-1, (s 2 - l)/(s 2 + 1)); (1 + xf' 2 (l - x)^ 2 ; N)} 



b 1= b 2 =N 
b 3 = 0, fa 4 = l/2 



b 1= b 2 =N 
b 3 =0, fa 4 = -l/2 



(3.3) 



where the final equality follows from (2.19). Recalling that Ziyi is defined as the RHS of (2.18), we 



see from the fact that Ziyi satisfies ( 2.17 ) that both cases of hyi in (3J3) satisfy the same differential 
equation. Comparing fl3.l| ) and ( |3.3D shows the r-functions in the latter also give the orthogonal and 
unitary symmetry gap probabilities, 

,2 



El (0; (-cj>, <j>);2N)=T~ hvi (-f TI ) 



(which is equivalent to a special case of (2.21)) and 



B 2 (Q;(-0,0);2iV) 



3.2 Bulk gap probability 



<C 2 + l) 



b 3 — 0, 64 = 1/2 



b 1= b 2 =N 
b 3 =0, fa 4 = l/2 



71 C 2 + l) 



b 1 =b 2 = N 

b 3 =Q, b 4 = -l/2 



(3.4) 



(3.5) 



In an obvious notation, the bulk scaled limit of (3.1) gives the formula 

(0;s) = -^ El (0,2s) + 



s) 



(3.6) 



Using the formula ( 2.22 ) for i^ ulk and a formula for E 2 ulk deduced from the analogue of ( |3.2| ), we have 
previously shown [Q0[ that this implies the Painleve transcendent evaluation 



E™*(0;s) = -^exp 



<J B (t) 



dt 

a=-l/2 i 



exp 



dt 

a=l/2 f 



(3.7) 



SO (7_b(<) 



a =-1/2 



and o b (t) 



where ob is specified by (2.22). Notice that the differential equation (2.23) is unchanged by a 1— > —a, 

differ in their characterization only by the boundary condition. The 



a=l/2 



definition ( [2.29 ) of T\w(t) and the characterization of any therein as the solution of ( 2.28| ) gives that 
(|3.7|) is equivalent to the r-function formula 



1 



£4 bulk (0;.s) = - WzV* 2 ) 



v 1 =v 2 = -l/2 



+ tuv(tt 2 s 2 ) 



v 1= v 2 = l/2 



(3.8) 



The boundary condition for any(t) when v\ = v 2 = — 1/2 is given by (2.31) while the corresponding 

sin 2\ft 



condition for V\ = V% = 1/2 is, from (2.24) 



(A ^ 
°ni' W ~ TT" 

t->o+ 2ir 



1 - 



2Vi 



t 3 / 2 



(3.9) 
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3.3 Soft edge scaling 



For the finite Gaussian ensemble the analogue of ( |3.l| ) is the coupled equations 12 



£; 2 (0;( S) oo);e- x ; 27V) 

= £a(0; (s, oo); e - ^ 72 ; 2N) \Ei(0; (s, oo); e - ^ 72 ; 2N + 1) + £a(l; (s, oo); e - ^/ 2 ; 2N + 1) 

+ £i(0; (s, oo); e^ 2 ; 2N + 1)^(1; (s, oo); e^ 2 ; 2iV) (3.10) 



£ 4 (0; (s, oo); e"^; iV) = J5 X (0; (s, oo); e^ 2 / 2 ; 2iV) + Ei(l; (s, oo); e"^/ 2 ; 2N + 1) (3.11) 

Here the only known quantity is Ei. In the soft edge scaling limit the number of distinct quantities in 
(|3.10|) is reduced, and one obtains the analogue of (|3.l|) |^2[ , 



(3.12) 



As noted in @, it follows from fl2.33|) and (|2.39|) that 

£f ft (0;( S ,oc)) - ^e^rO-^'W^e^r^^ + e-irt'-^WV^"*'*), (3.13) 



where q(t) satisfies ( 2.34 ) (this result was first derived in a direct calculation |3C[|). The first term in 
( |3.13| ) has in ( p. 57 ) been identified as a r- function. The second term differs from the first only in the 



sign of q{t). Since the differential equation (2.34) is unchanged by the replacement q t— > —q, we see that 



we can write the second term in (3.13) in a form formally identical to the first. Consequently 



E 



soft. 



(0;( S ,cx.)) = i{rW( S )+r( 2 )( S )} 



Q=0 



(3.14) 



where h in t^\s) is as in (2.57), while h in t^\s) is characterized as the solution of the same differential 



equation as in (2.57), but with the boundary condition 

h u (t) 



(iroo -iAiW + i{[Ai'(t)] 2 -i[Ai(i)] 2 } 



(3.15) 



which results by substituting (2.35) without the minus sign on the RHS in Q2.42D with a = 1/2 and 
recalling (|||) (c.f. (|1|)). 



3.4 Hard edge scaling 

In the case of the finite N Laguerre ensemble, the probabilities for the different symmetry classes of no 
eigenvalues in the interval (0, s) at the hard edge of the spectrum are related by coupled equations of the 
form ( 3.10| ), ( |3.11 ) fL2| , Consequently, in the scaled limit one obtains the analogue of (3.12) |Q 



£ 4 hard (0; (0, s )- a +!) = - + £^(0; (0, s); (a - l)/2) + 



hard / 



£ 2 hard (0;(0,s);a) 



and using (2.62) and the analogous result for E< 



hard 



£? ard (0;(0,s);(a-l)/2) 
we obtain @ 



(3.16) 



?(*) 



dt 



E l r d (0; (0, a); a + 1) = \ | cxp ( - ± jT ( log £ ) q 2 (t) dt) exp ( - ± jf ^_ 

+exp (^r( iog ^ 92(<) *) cxp (ir^ d< 



(3.17) 
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where q(t) satisfies (2.63). As with (3.13), the first term in (3.17) is the orthogonal ensemble result, which 
has been identified as a r-function in ( [2.77 ) above, while the second term differs from the first only in 
the sign of q(t). A further analogy is that ( 2.63 ), like ( 2.34 ) is unchanged by the mapping q i— > — q, so we 
have 



^(O; (0, s)i a + 1) = \ {t£ 3 (Vi) + rjg (^} 



r,=2 

Dl=-t)3 = -|(o-l) 

v 2 =— Vi=\{a+1) 



(3.18) 



where hy in t~ is as in (2.77), while hy in t£ satisfies the same equation (recall ( [2.63 )) but the 
boundary condition as determined by ( 2.65| ), ( 2.7C ), ( 2.76j ) and ( 2.64 ) is different because a minus sign 
must now be placed in front of (2.64). Explicitly we have for the second term 

xh v (x) ~ \xJ a {x) ~ -x 2 [J%(x) - J a+ i(x)J a -i{x)) (3.19) 

(c.f. ( 2.79| )). For each of the bulk, soft and hard edge cases, the r- functions in the evaluation of £4 give E\ 



and Ei according to formulas analogous to (3.4) and ( p.5| ). These are summarized in the accompanying 
table. 
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Gap Probabilities in Scaled Random Matrix Ensembles 



- Jardin d'Jimbo-Miwa-Okamoto 



Scaling Limit 


Orthogonal 


Unitary 


Symplectic 




Pup : vi = v 2 = ± l / 2 


P v : r? = 


= 1 


vi = v 2 = v 3 = v 4 = 


Bulk 


E^(0;(-s,s)) = t uv (tt 2 s 2 )\_ 1/2 


J B 2 bulk (0;(- S , S ))=T Iir (^V)|_ 
= T V (ins) 


.y 2 Tm'(7T 2 s 2 




£ 4 bulk (0; (-s/2, s/2)) = y 2 r Iir (^ 2 S 2 )|_i /2 + ^ruv^ 2 s 2 )\i/ 2 




ml ^ 

^iii'WL 1= , 2=T i/ 2 1 - 0+ ^ 


_^ sin2\/i 


t t 2 

Hy(t) ----- 






Pll 


a = 






Soft Edge 


E? 1 \0;(s,<x>))=t£\s) 


£r ft (0;( S ,oo))=T I ( I 1) (5)r I (2) (5) 

= -"(-2 1/3 «)L = _ 1/2 


sr(o;( s ,oo)) = y 2 4 1) ( s ) + y 2 r I (2) ( s ) 




4l' 2 \t) ~ ±y 2 Ai(t) + y 2 {[Ai'(i)] 2 - t[Ai(t)] 2 } -2^H n (-2^h)\ 2 ~ [Ai'(i)] 2 - t[Ai(t)] 2 




P v : V = 2 


vi = -v 3 = -y 4 (o - 1) u 2 = 


= -«4 = y 4 (a 
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Hard Edge 


E 1 hard (0;(0, S );(a-l)/2) = 4 1) (^) 


£ 2 hard (0; (0, s); a) = t v 1} (>/5) rf (^) 
= T„ r (s/4) 


£ 4 hard (0; (0, s); a + 1) = y 2 r v 1} {y/s) + ftr™ (v^) 




t/#' 2) (i) ^ TV 2 ^aW - V 4 i 2 [J 2 a(t) - Ja+l(t)Ja-l(t)] 
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